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/The  occurence  of  Intergranular  fracture  in  a  specific  situation  will 
v  depend  on  a  number  of  metallurgical  variables:  grain  size,  yield  strength, 
test  temperature,  and  notch  radius.  However,  in  all  reported  cases,  one 
observation  is  common.  Impurities  which  have  a  very  low  solubility  in  the  bulk 
have  segregated  to  the  grain  boundaries  and  weakened  them.  The  bulk  concen¬ 
tration  of  these  impurities  need  only  be  several  hundred  ppm  for  grain  boundary 
concentrations  of  5-10  atomic  percent  to  be  observed.  Examples  of  this  impurit 
induced  weakening  of  grain  boundaries  include  sulfur  in  iron  and  nickel, 
phosphorus,  tin,  and  antimony  in  steel,  bismuth  in  copper,  and  oxygen  in 
refractory  metals.  \ 
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i  If  the  environment  is  not  inert,. Intergranular  separation  may  occur  even 
though  the  fracture  mode  in  an  inert  environment  would  be  transgranular.  One 
much  studied  example  of  this  is  hydrogen  embrittlement. ^f  a  sample  is  tested 
either  in  hydrogen  gas  or  is  cathodically  charged  with  hytfrogen  and  then  testec 
in  an  inert  environment,  one  finds  that  the  mechanical  properties  of  the  material 
are  often  degraded  and  that  the  fracture  mode  has  become  more  intergranular. 

As-  with  impurity- induced  embrittlement,  many  additional  variables  can  affect 
the  degree  to  which  hydrogen  will  cause  embrittlement,  and  one  of  the  worst 
possible  cases  is  to  have  both  hydrogen  and  impurities  segregated  at  the  grain 
boundaries.  However,  it  is  clear  that  hydrogen  must  segregate  at  interfaces 
and  reduce  the  cohesive  strength  across  them. 

We  believe  that  an  Important  contribution  to  understanding  some  aspects 
of  these  complex  problems  of  embrittlement  in  metals  can  be  made  through 
theoretical  studies  of  bonding  in  model  systems.  The  philosophy  is  to 
isolate  groups  of  atoms  in  clusters  which  simulate  a  possible  local  environ¬ 
ment  of  interest  and  to  perform  detailed  quantum  mechanical  calculations 
using  molecular  orbital  theory  to  probe  the  nature  of  the  chemical  bonding. 
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INTRODUCTION  ! 

If  a  aetal  contains  very  few  impurities  and  Is  aechanlcally 
tested  In  an  inert  environment,  its  fracture  aode  will  be  either 
transgranular  cleavage  or  transgranular  ductile  alcrovoid  coales¬ 
cence.  If*  however,  the  aetal  Is  of  commercial  purity,  the  brittle 
fracture  aode  can  switch  froa  cleavage  to  Intergranular  separation. 
This  transition  Is  undesirable  because  intergranular  fracture  is 
generally  a  very  low  energy  process . 

\  ' 

The  occurrence  of  Intergranular  fracture  In  a  specific  situa¬ 
tion  will  depend  on  a  nuaber  of  metallurgical  variables:  grain 
sice,  yield  strength,  test  temperature,  and  notch  radius.  However, 
in  all  reported  cases,  one  observation  is  common.  Impurities  which 
have  a  very  low  solubility  in  the  bulk  have  segregated  to  the  grain 
boundaries  and  weakened  them.  The  bulk  concentration  of  these 
impurities  need  only  be  several  hundred  ppm  for  grain  boundary 
concentrations  of  5*10  atomic  percent  to  be  observed.  Examples  of 
this  impurity-induced  weakening  of  grain  boundaries  include  sulfur 
in  iron  and  nickel,  phosphorus,  tin,  and  antimony  in  steel,  bismuth 
in  copper,  and  oxygen  in  refractory  metals. 

If  the  environment  is  not  inert,  intergranular  separation  may 
occur  even  though  the  fracture  mode  in  an  inert  environment  would 


be  tr«asgr«ralir.  Ooe  much  studied  example  of  this  is  hydrogen 


__  embrittlement.  If  s  sample  is  tested  either  in  hydrogen  gas ,  or  is 

_  eathodically  charged  with  hydrogen  and  then  tested  in  an  inert 

environment,  one  finds  that  the  mechanical  properties  of  the  mate- 
'£  tial  are  often  degraded  and  that  the  fracture  node  has  become  more 
W  intergranular.  As  with  impurity-induced  embrittlement,  many  addi- 
_jL  tional  variables  can  affect  the  degree  to  which  hydrogen  will  cause 
embrittlement,  and  one  of  the  worst  possible  cases  is  to  have  both 
hydrogen  and  impurities  segregated  at  the  grain  boundaries.  Howevex; 
it  is  clear  that  hydrogen  must  segregate  at  interfaces  and  reduce  .  . 
_  .  the  cohesive  -strength  across  them  .  _ _ _ 

lie  believe  that  an  important  contribution  to  understanding 
some  aspects  of  these  complex  problems  of  embrittlement  in  metals 
.  .  can  be  made  through  theoretical  studies  of  bonding  in  model  systems. 
The  philosophy  is  to  isolate  groups  of  atoms  in  clusters  which 
simulate  a  possible  local  environment  of  Interest  and  to  perform 
--  detailed  quantum  mechanical  calculations  using  molecular  orbital 
theory  to  probe  the  nature  of  the  chemical  bonding. 

After  discussing  molecular  orbital -theory  and  the  Xa-scattered- 
wave  technique  in  the  next  section,  we  consider  in  the  following 
sections  two  recent  applications  of  this  approach  to  models  of 
potential  metallurgical  Interest.  * 


MOLECULAR  ORBITAL  THEORY  AND  THE  SCF-Xa-SW  CLUSTER  METHOD 

The  most  characteristic  feature  of  molecular  orbital  theory  is 
that  it  provides  an  independent-particle  interpretation  of  the 
electronic  structure  of  a  molecule  or  solid.  That  is,  each  electron 
of  the  system  has  associated  with  it  a  characteristic  spin-orbital 
and  energy  level  which  are  determined  by  solving  Schrodlnger's 
equation  for  an  electron  in  the  average  potential  energy  of  all  the 
other  electrons  in  the  system.  The  most  familiar  of  such  techniques 
is  the  Hartree-Fock  self-consistent  field  method,  in  which  each 
•  molecular  orbital  is  traditionally  represented  by  a  linear  combina¬ 
tion  of  atomic  orbitals  (the  HF-SCF-LCAO  method) .  For  all  but  the 
simplest  molecules  one  must,  at  least  for  the  present,  choose 
atomic-orbital  basis  sets  which  are  fairly  far  from  the  correct 
solutions  of  the  Hartree-Fock  Schrodinger  equations.  He  refer  to 
this  type  of  approximate  solution  of  the  Hartree-Fock  equations 
-  using  the  LCA0  representation  as  the  AHF-SCF-LCAO  method. 

-It  is  well  known  that  the  Hartree-Fock  theory  provides  a  poor 
description  of  the  electronic  structure  of  a  molecule  or  solid  at 
large  internuclear  distances  and  that  it  leads  to  an  Improper 
description  of  the  dlcsodatlon  of  such  systems.  In  order  to  cor¬ 
rect  for  this  behavior,  it  is  necessary  for  one  to  work  with  linear 
combinations  of  Slater  determinants  of  molecular  orbitals,  the 


3 


so-called  "configuration- interaction"  (Cl)  procedure,  which  is  not 
only  conduce  ted  computationally ,  hut  also  destroys  the  independent¬ 
ly-  olectron  interpretation.  4  further  difficulty  is  that  in  performing 
numerical  calculations  using  the  LCAO  approximation  at  either  the 
•  r  4BF  or  Cl  levels  of  sophistication,  the  nunber  of  two-electron 
■,  -  multicentcr  integrals  which  must  be  calculated  even  for  moderately 

sized  systems  becomes  enormous.  For  example,  in  a  recent  calcula- 
don  on  a  molecule  consisdng  of  26  atoms,  all  of  atomic  numbers 
^ .  less  than  ten,  10 ^  Integrals  had  to  be  evaluated  and  the  calcula- 
...  dons  required  the  equivalent  of  192  h  on  an  IBM  360/195  computer, 

_ pne  of  the  fastest  machines  available.  _ ! _ 


*  A  much  more  practical  independent-electron  theory,  which  avoids 

~  most  of  the  difficulties  associated  with  the  4HF-SCF-LCA0  and  Cl 
^  .  methods,  is  the  recently  developed  SCF-Xn-SW  cluster  method.  ~Based  .  ✓ 

on  the  combined  use  of  Slater *s  Xa  statistical  theory  of  exchange 
:  correlation  and  Johnson's  multiple-scattered-wave  formalism,  this 
technique  allows  one  to  calculate  accurate  self-consistent  spin- 
,1.:  orbital  energies ,  wave  functions ,  charge  densities ,  and  total 
energies  for  complex  molecules  and  clusters,  yet  requires  only 
relatively  small  amounts  of  computer  time.  Like  the  AHF-SCF-LCAO 
approach,  the  SCF-Xn-SW  method  is  an  approximate  theory  of  many- 
electron  systems,  but  It  Is  not  an  approximation  to  Hartree-Fock 
theory.  First  of  all,  in  contrast  to  HF  theory,  the  Xa  statistical 
expression  for  the  total  many-electron  energy  automatically  goes 
•  to  the  proper  separated-atom  limit  as  the  lnternudear  distances 
are  increased  to  infinity,  l.e.,  as  the  system  dissociates.  Second, 

It  has  been  shown  that  the  Xa  theory  satisfies  Fermi  statistics, 

: thereby  ensuring  the  proper  ordering  and  occupation  of  electronic 
energy  levels,  while  the  HF  theory  does  not.  Third,  it  has  been 
proven  that,  under  the  proper  conditions,  the  Xa  statistical  total 
energy  rigorously  satisfies  both  the  virlal  and  Hellmann-Feynman 
theorems,  which  facilitates  the  calculation  of  the  equilibrium 
cohesive  energies  of  molecules  and  solids.  Finally,  the  SCF-Xa-SW 
method,  in  conjunction  with  the  "transition-state"  concept,  leads  ^ 
to  an  accurate  approximate  description  of  the  excited  electronic 
states  and  charge  distributions  of  polyatomic  systems,  including 
the  effects  of  orbital  relaxation  which  are  usually  neglected  in 
applying  Koopmans  theorem  in  the  AHF-SCF-LCAO  approach. 

The  SCF-Xa-SW  method  and  transition-state  procedure  have 
already  been  applied  successfully  to  a  vide  range  of  molecules  and 
clusters  in  solids.  Because  the  required  computational  effort  does 
not  increase  inordinately  with  the  number  of  electrons  per  atom  or 
with  the  nuober  of  atoms  in  the  molecule  or  cluster,  the  SCF-Xa-SW 
method  is  ideally  suited  for  describing  the  electronic  structures 
of  transition-metal  complexes  and  clusters  representing  the  local 
molecular  environment  of  a  grain  boundary. 


.WLBCULAR  OUITAL  MODEL  IDE  GEAIH  BOUNDARY  EMBRITTLEMENT  .OF .NICKEL 
SYSULFUR 

As  discussed  2a  tit*  last  section,  the  theoretical  aethod  we 
■he  to  investigate  the  electronic  structure  of  the  local  environ- 
■ent  at  a  grain  boundary  constats  of  two  parts .  The  first  is  to 
.use  a  cluster  of  atone  to  represent  the  local  environment  end  the 
second  is  to  use  molecular  orbital  theory  to  solve  for. the  elec-.  . 
tronic  structure  of  the  cluster  of  atone.  _ _ _ 


To  apply  this  technique  to  grain  boundaries  one  wist  first  . 

choose  a  cluster  that  reseobles  in  some  way  the  local  environment 

an  impurity  aton  might  have  were  it  at  the  grain  boundary.  The _ 

choice  of  such  dusters  is  considerably  simplified  due  to  the  .  .. 

recent  work  of  Ashby,  Spaepen,  and  Williams,  who  have  shown  that 
the  structure  of- grain  boundaries  la  fee  solids  can  be  described  in 
terms  of  atoms  at  the  vertices  of  Bernal  deltahedra.  The  individual 
deltahedra  are  small  clusters  containing  four  to  thirteen  atoms  and 
thus  provide  a  natural  initial  model  of  local  atomic  structure  with 
which  to  investigate  the  electronic  properties  by  SCF-Xa-SW  calcu¬ 
lations.  In  order  to  investigate  the  feasibility  of  applying  such 
methods  to  this  problem  we  have  chosen  to  use  one  particular  cluster 
in  this  work. 

The  particular  duster  we  Chose  Is  the  tetragonal  dodecahedron 
shown  in  Fig.  1.  It  contains  eight  atoms  and  has  D2d  point  group 
symmetry.  We  do  not  claim  that  this  structure  is  necessarily  a 
prevalent  one  in  grain  boundaries.  However,  it  does  contain  four 
atoms  (nusbered  1,  2,  3,  and  4)  slightly  but  syonetrically  displaced 
from  the  x-y  plane  and  four  other  atoms  further  displaced  from  this 
•plane.  One  can  roughly  consider  that  the  four  atoms  nearest  the 
x-y  plane  are  grain  boundary  atoms  and  that  the  other  four  are  in 
the  first  atomic  layers  away  from  the  precise  boundary. 

i  It  is  well  established  experimentally  that  sulfur  enfcrlttles 
nickel.  In  order  to  investigate  the  electronic  effects  of  a  sulfur 
atom  at  a  grain  boundary  in  Hi,  we  have  carried  out  two  calculations 
on  a  model  system  in  which  the  Hi  atoms  were  placed  at  each  of  the 
eight  sites  of  the  gluster  shown  in  Fig.  1  with  an  internudear 
separation  of  2*49  A.  In  the  first  calculation,  we  considered  the 
cluster  containing  only  the  eight  nickel  atoms .  For  the  second 
calculation,  a  sulfur  atom  was  placed  at  the  center  of  the  cluster. 
An  interstitial  site  of  this  type  is  a  possible  position  Which  a 
sulfur  atom  would  choose  at  a  grain  boundary.  For  these  two  con¬ 
figurations  we  have  determined  the  molecular  orbitals  for  the 
clusters,  their  corresponding  energy  levels,  and  the  valence  charge 
densities. 

-  •  The  energy  levels  for  teach  cluster  are  shown  in  Fig.  2.  We 
note  that  in  the  MlgS  cluster  several  of  the  energy  levels  have 
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Figure  2 


been  lowered  subs  tea  daily.  .  Those  Chet  have  been  lowered  have  . 
orbital  wave  functions  which  are  highly  localized  on  the  sulfur 
atoa  and  its  four  nearest  neighbors,  i.e.,  the  four  nickel  atoas 
closest  to  the  x-y  plane.  Based  on  this  slaple  analysis  one  night 
assune  chat  sons  Ni-Ni  bonds  have  been  weakened.  n‘ 

f  'a  *•  * 

The  bonding  nay  be  better  understood  if  one  analyzes  the 
aolecular  orbitals  of  the  clusters  in  more  detail;  however,  in  .  _ 
this  particular  case  a  aore  succinct  picture  of  the  bonding  in  the 
cluster  can  be  obtained  by  sunning  the  orbital  charge  densities 
over  all  the ‘occupied  valence  orbitals.  One  can  then  plot  contours 
of  the  total  valence  charge  density  which  exists  on  any  plane. 

These  are  shown  in  Fig.  3  for  the  x-z  plane.  In  this  figure  we 
have  plotted  the  sane  set  of  contours  for  both  clusters.  The  three 
outermost  contours  (nunbered  k,  2,  end  3)  in  both  clusters  surround 

all  atoas.  However,  contour  number  4  (which  is  shown  with  a  _ 

heavier  line  for  clarity)  is  quite  different  in  the  two  clusters. 

In  Nig,  it  shows  bonding  between  the  pairs  of  nickel  atoms  1  and  7, 

7  and  8,  and  4  and  8,  where  the  characteristic  "neck”  regions 
between  the  pairs  of  atone  is  responsible  for  bonding.  In  NigS, 
the  shape  of  contour  number  4  is  quite  different.  The  character¬ 
istic  neck  bonding  region  between  atone  1  and  7,  and  4  and  8  is 
gone  and  consequently  so  is  the  concomitant  strong  bonding  between 
the  two  types  of  nickel  atoms.  Furthermore,  the  size  of  the  neck 
region  between  atoms  7  and  8  is  reduced  from  that  in  Nig,  implying 
a  weakening  of  this  metal-metal  bond  as  well.  Contour  4  in  NigS 
is  split  into  two  distinct  regions.  The  upper  region  of  contour  4 
and  an  inner  contour  of  the  nickel  atoms  nearest  the  sulfur  atom 
enclose  both-  the  nickel  atoms  and  the  sulfur.  This  would  lead  to 
very  strong  bonds  between ‘tills  set  of  nickel  atoms  and  the  sulfur. 
However,  there  is  no  such  bonding  interaction  between  the  sulfur 
atoms  and  the  other  nickel  atoms. 

From  these  results  a  clear  picture  begins  to  emerge.  Adding 
sulfur  to  the  cluster  of  nickel  atoms  causes  strong  NiS  bonds  to 
be  formed  between  the  nickel  atoms  nearest  the  sulfur.  Yet  this 
weakens  bonds  which  the  nickel  atoms  would  ordinarily  have  with 
the  other  nickel  atoms  that  are  not  as  dose  to  the  sulfur  atom. 

In  tills  cluster  the  strong  Ni-S  bonds  involve  the  sulfur  atom  and 
the  four  nickel  atoms  nearest  the  x-y  plane.  We  suggested  earlier 
that  these  atoms  might  be  thought  of  as  being  in  the  grain  boundary. 
However,  the  bonds  between  the  nickel  atoms  nearest  the  x-y  plane 
and  those  farther  away  are  weakened.  These  bonds  might  represent 
those  from  the  boundary  out  to  the  first  layer  of  the  bulk.  Clearly 
if  those  are  weakened,  the  stress  required  for  fracture  would  be 
less. 

Therefore,  based  on  the  calculations  we  can  suggest  an  elec¬ 
tronic  mechanism  by  which  grain  boundary  embrittlement  could  occur. 
We  stress  that  these  are  the  results  from  calculations  on  only  one 


system.  Furthermore,  we  have  chosen  sa  embrittling  element  which 
probably  occupies  an  interstitial  site  at  the  grain  boundary. 

Other  eabrittiing  eleaents  would  clearly  be  substitutional.  Also, 
•any  other  types  of  geometrical  configuration  are  possible  at  the 
boundary.  Yet,  with  these  results,  we  believe  that  it  la  not  pre¬ 
sumptuous  to  propose  that  the  above-described  mechanism  is  one 
mechanism  by  which  segregated  elements  embrittle  the  grain  boundary. 

MOLECULAR  ORBITAL  MODEL  FOR  HYDROGEN  EMBRlITLEKENT  OFIRON 

The  causes  and  mechanisms  of  hydrogen  embrittlement  have  been 
and  will  continue  to  be  the  subjects  of  detailed  Investigation. 

As  a  result,  several  models  have  been  offered  in  explanation  of  the 
embrittlement  .phenomenon.  These  models,  however,  cannot  account  . 
for  embrittlement  at  the  atomic  level,  as  a  chemical  interaction 
between  hydrogen  and  the  embrittled  metal.  Zt  is  the  purpose  of 
this  paper  to  present  die  results  of  SCF-Xa-SW  molecular-orbital 
calculations  on  iron-hydrogen  clusters  in  order  to  explain  hydrogen 
embrittlement  of  ferrous  alloys  at  the  most  fundamental  level. 

Hydrogen  atoms  are  most  frequently  found  in  octahedral  holes 
of  transition  metals.  If  molecular  hydrogen  were  adsorbed  within 
tetrahedral  holes  then  concerted  dissociation  of  molecular  hydrogen 
would  result  in  placing  atomic  hydrogen  in  octahedral  holes. 
Accordingly,  the  four-atom  iron  cluster  shown  in  Fig.  4  will  be  used 
to  model  the  tetrahedral  hole  and  the  site  for  hydrogen  adsorption 
in  bcc  iron.  (Both  the  tetrahedral  and  octahedral  holes  of  bcc 
iron  are  misnamed  since  they  actually  present  a  D2<1  and  tetragonal 
symmetry  respectively.)  In  an  earlier  paper,  the  results  of  SCF- 
Xa-SW  molecular-orbital  calculations  of  four-,  nine-,  and  fifteen- 
atom  clusters  of  iron  were  presented.  It  was  shown  that  these 
clusters  provide  a  model  for  the  local  electronic  structure  of  bcc 
iron.  A  comparison  of  these  results  indicates  that  both  the  four- 
atom  and  fifteen-atom  clusters  exhibit  the  large  exchange  splitting 
and  high  magnetic  moment  characteristic  of  bulk  iron.  Net  spin 
electron  nunber  per  atom  for  Fei5  is  2.7  and  for  the  Fe4  cluster 
2.5  compared  with  2.12  for  the  bulk.  The  density  of  states  in  the 
four-atom  cluster  is  sufficiently  similar  to  the  density  of  states 
in  the  15-atom  cluster  to  justify  the  use  of  Fe4  as  a  model  for  a 
site  of  hydrogen  adsorption.  It  is  important  to  note  that  the 
fifteen-atom  cluster  provides  three  environments  for  the  atoms  of 
the  cluster.  In  bulk  iron  (neglecting  surface  effects)  each  iron 
atom  Is  in  the  same  environment.  It  is  this  fact  which  is  respon¬ 
sible  for  the  veriation  of  properties  predicted  by  the  cluster 
method  with  those  measured  in  bulk  iron.  In  a  calculation  of  this 
type,  however,  where  an  impurity  is  being  modeled,  the  impurity 
breaks  the  symmetry  of  the  bulk  and  there  are  a  nunber  of  environ¬ 
ments  around  the  impurity  site.  As  a  result,  the  cluster  method  is 
ideally  suited  for  the  study  of  dilute  impurity  states  and 
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consequently  smaller  clusters  then  those. used  In  the  study  of  bulk 
properties  can  be  used.  -  .  -  .  _L 

The  eolecular-orbi  tel  energy  levels  of  s  hydrogen  aolecule 
Inserted  into  the  four-atom  cluster  is  shown  In  Pig.  S.  The  a^ 
and  b2  orbitals  of  the  four-atom  Iron  cluster  have  the  a  erne  sym¬ 
metry  as  the  hydrogen  aolecule  bonding  and  antibonding  orbitals. 
.These  hydrogen  aolecular  orbitals  couple  with  the  lowest  lying  j. 
ai  and  b2  orbitals  of  the  Pe4  cluster  to  produce  the  aolecular 
oxbitals  shown  In  Figs.  6  .and  7.  As  a  result  of  this  coupling  the 
lowest  lying  ai  and  b2  orbltala  of  Fe4H2  are  stabilized  relative, 
to  the  position  of  these  orbitals  In  Fe4.  The  a^  orbital  Is  pre¬ 
dominantly  a  B-B  bonding  orbital  the  B-Fe  interaction  is  non¬ 
bonding.  The  b2  orbital  In  the  Fe4  cluster  is  d-d  it -bonding  between 
nearest  neighbors  and  6-bonding  between  second-nearest  neighbors. 

The  addition  of  B2  strengthens  the  bonding  Interaction  between 
second-nearest  neighbors,  has  little  effect  on  the  nearest-neighbor 
interaction,  and  is  strongly  antibonding  between  the  hydrogen  atone . 
The  B-B  distance  in  this  cluster  is  the  molecular  B-B  distance. 

This  is  the  distance  which  is  realized  when  there  is  no  electron 
density  in  a  B-B  antibonding  orbital.  In  this  cluster  some  of  the 
electron  density  from  the  metal  cluster  b2  orbital  is  transferred 
into  the  antibonding  orbital  of  the  hydrogen  molecule.  The  hydrogen 
molecule  fl-B  distance  is  no  longer  a  stable  equilibrium  distance 
and  the  hydrogen  atoms  will  begin  to  dissociate. ■  All  transition 
metals  are  known  to  dissociate  hydrogen  molecules.  The  efficiency 
of  this  dissociation  will  depend  on  the  amount  of  overlap  between 
the  metal  orbital  of  suitable  symmetry  and  the  antibonding  orbital 
of  the  hydrogen  molecule.  Other  than  the  low-lying  a^  and  b2 
orbitals,  no  other  orbitals  of  this  cluster  have  any  significant 
hydrogen  admixture.  Two  additional  spin  orbitals  are  occupied, 
however,  by  the  two  electrons  contributed  by  the  hydrogen  molecule. 

Figure  8  gives  the  molecular  orbital  energy  levels  for  the 
same  arrangement  of  iron  atoms  but  now  the  H-H  distance  has  been 
Increased  as  a  consequence  of  the  electron  density  In  the  H-H  . 
antibonding  orbital.  Again  the  lowest  lying  a^  orbital  Is  the  only 
orbital  with  significant  H-H  bonding  contributions  (Fig.  9).  The 
lowest  lying  b2  orbital  is  H-H  antibonding,  as  in  the  previous 
case,  but  dlls  Interaction  has  been  lessened  as  is  clearly  seen  in 
Fig.  10.  There  is  a  reasonably  strong  H-Fe  bonding  interaction 
within  dlls  b2  orbital,  and  this  interaction  should  be  maximized 
when  the  hydrogen  atom  is  sitting  between  the  Iron  second-nearest 
neighbors,  l.e.,  when  the  hydrogen  atoms  are  In  the  tetragonal 
holes  of  the  bcc  structure. 

While  hydrogen  adsorption  and  dissociation  may  be  explained 
in  a  fairly  straightforward  manner  by  molecular-orbital  theory, 
attempts  to  explain  embrittlement  at  the  molecular  level  must  be 
assisted  by  considering  embrittlement  as  a  competitive  process 


between  two  aechaaisas.  The  fleet  of  these  Is  the  process  of 
piss  tic  flow  facilitated  by  dislocation  aoveaent.  The  second  is 
characterised  by  nuclear  notion  in  a  direction  which  will  counter 
die  strain  energy  of  the  crystal,  a  purely  elastic  response-.  An 
energy  of  activation  nay  be  associated  with  each  of  these  processes . _ 
The  activation  energy  for  the  first  process  will  depend  predomi-  _ 
nately  upon  the  structure  of  the  crystal,  the  Burgers  vector,  etc., 
while  the  activation  energy  of  the  second  process  is  slnply  the  ..  — 
fracture  energy  of  the  crystal.  A  ductile  netal  responds  to  .  ... 
excessive  strain  through  the  first  aechanlsn,  and  this  nay  be 
attributed  to  the  fact  that  the  activation  energy  for  plastic  flow 
la  significantly  less  than  the  activation  energy  for  brittle  frac¬ 
ture.  A  brittle  naterial  will  respond  to  strain  by  fracturing 
before  plastic  ~flow  begins  because  the  activation  energy  for  dis¬ 
location  aoveaent  is  higher  than  the  activation  energy  for  brittle 
fracture.  A  ductile-brittle  transition  aust  be  the  result  of  a 
aechanlsa  that  either  raises  the  activation  energy  for  plastic  flow 
or  lowers  the  energy  for  brittle  fracture.  — - -  .... 

II  the  process  of  hydrogen  adsorption  affects  the  activation 
energy  for  dislocation  flow,  then  hydrogen  adsorption  aust  alter 
the  structure  of  bcc  iron.  Structural  changes  are  often  aanifest 
as  changes  in  the  aolecular  orbital  energy  distribution  near  the 
Feral  level.  Figure  11  shows  the  aolecular  orbital  energy  distribu¬ 
tion  near  the  Feral  level  for  a  tetrahedron  of  iron  atoas,  the 
four-atoa  cluster  aodeling  bcc  iron,  and  three-hydrogenated-lron 
cluster  corresponding  to  adsorbed  aolecular  hydrogen,  adsorbed 
dissociated  aolecular  hydrogen  and  a  hydrogen  saturated  cluster, 
Fe4H£.  In  the  tetrahedral  cluster,  there  is  an  unoccupied  t2 
orbital  near  the  Ferad.  level,  this  orbital  is  split  when  the  syw- 
aetry  of  the  cluster  is  reduced  froa  tetrahedral  to  D2d  as  in  the 
4-aton  bcc  cluster.  The  e  orbital  is  stabilised  and  becoaes  occu¬ 
pied  while  the  b2  orbital  is  destabilised.  This  process  is  real- 
niscent  of  the  Jahn-Teller  effect  of  inorganic  coaplexes,  where  a  ^ 
coaplex  spontaneously  distorts  in  order  to  lower  its  total  energy. 
This  leads  to  the  conjecture  that  the  bcc  structure  of  iron  is  a 
consequence  of  a  Jahn-Teller  distortion  of  a  fee  structure.  When 
hydrogen  aolecules  are  added  to  the  bcc  cluster  the  e  orbital  is 
depopulated  and  is  nearly  degenerate  with  the  b2  orbital,  and  is 
slallar  to  the  energy  distribution  of  the  tetrahedral  iron  cluster. 
The  dissociation  of  the  hydrogen  aolecule  stabilises  the  b2  orbital 
but  the  addition  of  aore  hydrogen  seeas  once  again  to  depopulate 
this  orbital.  Hydrogen  adsorption  would  appear  to  proaote  a  dose- 
packed  phase  transition  near  the  site  of  adsorption.  There  is 
evidence  that  hydrogen  adsorption  does  proaote  phase  transforaa- 
tlons .  In  HCP  Wtsab,  hydrogen  adsorption  is  known  to  lower  the 
stacking  fault  energy  indicating  that  hydrogen  is  proaoting  the 
fcc-hcp  phase  transition. 

This  la  Just  the  type  of  effeet  that  would  raise  the 


~~  actlTatloa  n«Z0  for  dislocation  aovtacat  aroiod  th«  sites  of  *  * 

hydrogen  adsorption.  However, 'no  surface  reconstruction  of  Aron 
due  to  hydrogen  adsorption  has  ever  been  noted  and  fee  aetala, 

^notably  nickle,  eabrlttle  In  a  hydrogen  environment .  While  hydrogen 
undoubtedly  has  an  effect  on  die  structure  of  bee  Iron,  It  la  unlike- 
ly  that  this  alone  can  account  for  hydrogen  eabrlttlenent.  ‘ 

._U  ...  .V  » 

'Jt  ...  .  This  leaves  the  possibility  that  hydrogen  adsorption  lowers _ 

die  activation  energy  of  brittle  fracture.  If  this  Is  the  case 
then  hydrogen-adsorbed  clusters  should  have  aolecular  orbitals  that 
show  siallarities  to  aetal  clusters  that  have  been  strained,  for 
such  clusters  require  less  additional  energy  to  fracture,  i.e.,  the 
activation  energy  for  brittle  fracture  has  been  lowered.  Accordingly, 
the  aolecular  orbital  energy  levels  for  the  four-aton  bcc  cluster 
in  which  the  Iron-  atoas  have  been  strained  20Z  along  the  twofold 
rotation  axis  was  calculated.  The  relative  energy  distribution  of 
the  najority  spin  states  for  this  cluster  la  shown  In  Fig.  12, 
along  with  the  four-aton  nodal  of  bcc  iron  and  the  dissociated  ...  . 
'aolecular  hydrogen-iron  clusters.  (The  ainorlty  spin  states  for 
all  clusters  show  the  sane  ordering  of  the  energy  levels  as  do  the 
najority  spin  states.  Ho  inforaation  Is  lost  by  studying  only  the 
najority  spin  states,  while  the  forest  of  levels  is  considerably 
thinned.) 

The  process  of  straining  the  bulk  naterlal  results  In  con¬ 
siderable  rearrangenent  of  the  energy  levels  within  the  energy 
aanifold.  A  rough  neasure  of  the  difference  between  two  energy 
nanifolds  is  the  nuaber  of  Intersections  between  the  connecting 
lines  of  the  nanifolds.  There  are  ten  intersections  between  the 
bulk  and  the  strained  naterlal,  while  there  are  only  two  intersec¬ 
tions  connecting  the  strained  and  the  hydrogenated-iron  energy 
nanifolds.  Saturation  of  the  Fe4  cluster  with  hydrogen  results  in 
alnost  no  rearrangenent  of  the  energy  nanifold.  What  rearrangenent 
occurs  is  confined  to  the  lowest  five  orbitals  of  the  unifold. 

"s 

The  direction  in  Which  the  nolecular-orbltal  energy  levels  move 
upon  straining  can  be  seen  to  be  the  result  of  removing  electron 
density  from  between  nearest  neighbors,  and  replacing  it  in  the 
region  perpendicular  to  the  direction  of  strain,  i.e.,  between 
second-nearest  neighbors..  The  result  should  be  a  contraction  of 
the  cluster  perpendicular  to  the  direction  of  strain.  This  is  the 

effect.  Within  this  analysis  it  is  clear  that  hydrogen  atoms 
7  in  tetragonal  holes  will  mimic  the  behavior  of  strain  by  increasing 
7  the  electron  density  between  second-nearest  neighbors.  It  is  worth 
noting  that  only  the  two  lowest  lying  orbitals  of  the  hydrogenated- 
iron  cluster  have  any  hydrogen  admixture;  therefore,  the  rearrange- 
unt  of  the  other  orbitals  must  be  the  result  of  a  relaxation  effect 
.due  to  the  stabilisation  of  the  lowest  a^  and  b2  orbitals. 

It  seems  reasonable  then  that  if  the  lowest  lying  aj  and  b2 


